Equation of state for hard sphere fluids with and without Kac tails 
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In this note, we propose a simple derivation of the one dimensional hard rod equation of state, 
with and without a Kac tail (appended long range and weak potential). The case of hard spheres 
in higher dimension is also addressed and it is shown there that our arguments -which avoid any 
mathematical complication- allow to recover the virial form of the equation of state in a direct way. 

I. INTRODUCTION 

The equation of state is a central object in the statistical mechanics description of matter. Bearing the signature 
of inter-particle interactions and being easily measured experimentally, it establishes a connection between micro- 
scopic and macroscipic features that long remained elusive. [1] In particular, many efforts in the 19th century aimed 
at explaining the deviations from the ideal gas law. Progresses in the field have been hindered by the positivist 
storm of criticisms directed against kinetic theory and the associated atomistic viewpoint. Today, the machinery of 
statistical mechanics offers a consistent framework for those problems, |2i, iS] with the pedagogical drawback that very 
few interacting systems have an equation of state that can be written in closed form. 

The one dimensional hard-rod fluid (also called Tonks or Jepsen gas [1, @) is one notable exception. For N rods 
with a distribution of lengths {^i}i<i<N ^ enclosed in a line of total length L, the pressure P can simply be written [3| 

P - f^, (1) 
1 - 77 

where kT is the thermal energy, p — N/L denotes the density, and rj = J^i^i/^ ^^e line covering fraction (so that 
T]/ pis the mean rod length). In addition, several equilibrium and non-equilibrium quantities can be computed, [sj and 
the model provides a reference system for perturbative treatments, one of which being addressed below, with inclusion 
of a so-called Kac pair potential of interaction between the particle s, [g, H 

In standard textbooks [1, Q as well as in the original papers. P. [sl. lol. Il0| the derivation of Eq. ([1]) is not straightfor- 
ward, and it is our purpose here to propose an alternative concise argument relying on simple physical considerations 
rather than partition functions and mathematical computations. When transposed to systems of higher dimensions 
d > 1, the argument provides the exact equation of state of a hard sphere fluid, which however cannot be written in 
closed form as an explicit function of density. [lH 

II. EQUATION OF STATE OF HARD ROD AND HARD SPHERE FLUIDS 

We begin with the one dimensional case (hard rods). To achieve our goal of computing exactly the pressure, we 
consider several related "Gedankenexperiments" . We first make the remark that upon rescaling all lengths in the 
problem by a factor 1 + e 

+ e% for all 1 < i < A^, L ^ (1 + e)i, (2) 

the reversible work (5Wtotai associated with this process reduces to its ideal gas contribution since excluded volume 
is irrelevant in this transformation (it is of course essential to rescale the "box" length L) ; only the ideal entropy of 
mixing is affected in this simple "zoom" . We then have 

^Wtotai = -pkT6L, (3) 

where 5L = eL is the total length change. Equation ^ holds beyond hard body interactions and applies whenever the 
problem at hand is governed by a unique length scale, [l^ It would in particular hold for Lennard- Jones interactions. 
For 5L > 0, the work ^Wtotai is negative. If only the box length is expanded, the work (5Wtotai received by the gas is 
obviously negative. The additional particle expansion considered here amounts to a positive contribution to (5Wtotah 
which turns out to be smaller in absolute value than the former one. The resulting cost is negative, as the following 
discussion explicitly shows. 
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We now perform the expansion ([2]) in two steps and compute separately the corresponding reversible work required 

a) First, rod sizes are slowly and sequentially rescalcd (i.e. one at a time): 

£i^(l + e)^i, ^2^(l + e)^2... 

b) Second, the confining box size is expanded L — > (1 + e)L. 

It is essential to realize that in step a), any rod that is expanded behaves as a confining wall (piston) which "pushes" 
the fluid so that the work needed to rescale particle i is P5ti — Peii (we comment further below this fundamental 
fact). Consequently, the work needed to perform step a) reads 

N 

SWa = J2PS£^ = PSLp, (4) 
1=1 

where Lp is the total length occupied by the rods (here, SLp ~ rj6L). For step b), the work reduces to 

SWb = -P6L. (5) 
Summing both contributions, we have JWtotai = SWa + SWb and gathering results, this implies 

- pkTSL = Pr]SL- P6L =^ P = (6) 

1-77 

which is the correct result. We stress that irrespective of the sign of the scaling parameter e, the moves considered in 
steps a) and b) do not produce overlaps between particles nor between particles and confining boundaries, in the same 
way as a moving boundary like a piston does not lead to overlaps. We also note that from the form of the equation 
of state dl]), all virial coefficients are unity. This can be checked by explicit calculations, 

Following the same line of reasoning, we investigate the problem of hard spheres in dimension d > 1. One still has 

SWtot.1 ^ SWa + SWb with SWtot.i = ~pkTSV, (7) 

V being now the volume of the system, and SWb — —PSV. However, SWa is not as straightforwardly related to the 
pressure as in ^ and to obtain its expression, we introduce the pair correlation function ofr).[3l. For the sake 
of simplicity, we restrict to the monodisperse case where all particles have the same diameter a. The force per unit 
area felt by a given particle is of kinetic origin and reads pkTg{a) ; remembering that a particle is surrounded by an 
excluded volume sphere of diameter 2a where no other particle's center of mass may be found, we can write the work 
needed to expand a into a -\- Sa as 

SWa^a+Scr = pkT g{a) (5V;wcep, (8) 

where i5V^wocp is the volume change of the excluded volume sphere. The latter quantity is related to the d-dimensional 
surface Sdi'ia) of a sphere with radius 2cr through 

'JV^swccp = Sa{2a) ^ = 2'^-i5rf(a) ^. (9) 
Summing over all particles in the system, we get 

SWa = SW^^,+sa = pkTg{a)2''-^ SVp (10) 

i 

where, keeping previous notatioii, SV^j is the change of the total volume occupied by the spheres and therefore 
SVp/SV — 77, the volume fraction. [1J| Returning to Eq.(l7]), we obtain 

-pkT = pkT2''-^Tig{cj)-P, (11) 

and hence the equation of state 

= l + 2'^-Sg(a). (12) 

To our knowledge, the simplest derivation of this result involves the virial theorem [lT| and turns out to be more 
complicated and physically less transparent. [l^ As a by-product, comparing (|12p and ^ in one dimension, we obtain 
that the pair correlation function at contact for hard rods takes the value 

g(^) = 1/(1 - 77), (13) 
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a standard result. We also emphasize that it is straightforward to generalize (|10p and ((T2l) to the polydisperse 
case with size distribution /(cr) (see e.g appendix B of Ref. [IBl), with the result 

^ -l=rj fdad<j'f{a)f{a')g{a/2 + a'/2)^^^±^^, (14) 



pkT 'J ' ' ' (ct-^) 

where (ct'') / a'^f{a)da. 

III. INCLUSION OF A KAC TAIL 

It proves instructive to consider also the case of particles interacting with a so-called "Kac tail" since not only can 
the proper pressure be recovered, but also some light be shed on the nature of the expansion processes underlying our 
arguments. We assume that in addition to the usual hard sphere term, particles interact with a very long range and 
weak pair potential so that the interaction potential reads 

d(r) = for r<(T = £ , . 

^ ' \ 7exp(— r/r*) for r > a, ^ ^ 

where the range r* is larger than any microscopic distance {pr* » 1) and the amplitude 7 is small. [13] This quantity 
can be positive (repulsive tail) or negative (attractive tail). It has been shown -with all mathematical is dotted- 
that the corresponding equation of state is of the van der Waals form and reads Q 

P=^ + ap\ (16) 

1—7/ 

where a = 7 r*. This result has been generalized to an arbitrary space dimension 7] in which case again, the correction 
to the hard sphere pressure is simply ap^, with now ct = J (/"(r) d'^r/2, where the integral runs outside the excluded 
volume sphere (i.e. r > a). Actually, the result holds irrespective of the precise form of the interaction potential, Q 
as is also clear from the following argument. We restrict to c? = 1, but space dimension is largely immaterial. From 
SWtotai = SWa + 5Wb, we have 

{P - pkT)SL ^ SWa, (17) 

and our objective is now to compute the latter term. In step a), it is understood that the range of the potential r* 
is expanded, at fixed amplitude 7, to remain proportional to the rod size £: 6r* /r* = 6£/£ = 6L/L. A key point in 
the analysis is that the tail r > i oi potential does not affect the relative positions of the particle from what they 
would have as hard rods at the same density p.Q This follows from the fact that the amplitude 7 should be taken 
extremely small. [TTj The pair distribution function g{r) is thus unaffected by the tail and we have in particular the 
contact value (jlSp. When a given rod of size £ is expanded, the work necessary can be written as the sum of a kinetic 
contribution (one has to "push" the particles that are in direct contact with the particle of interest), and a (long 
range) contribution arising from the tail of the potential : 

SWe-.e+Si = Pcont.ctkT5£ + SWt.n - S£ + SWt.n. (18) 

1 - T] 



Since any particle experiences an average potential energy 2^r*p, as follows from integrating (|15p . we have (SVFtaii = 
6(2jr* p) where the variation is computed at fixed amplitude and density. Summing over all particles, we obtain 

SWa =^vSL + ^N2jpSr*, (19) 

where the factor 1/2 corrects for double-counting. Remembering that 6r* /r* — 6L/L, we arrive at 

SWa - Y~ VSL + yr*^ p^ SL. (20) 



Inserting this result into P?|) yields Eq. a result otherwise difficult to derive. 0] As is intuitively clear. 

Equation (fTB|) indicates that repulsive interactions (7 and a positive) lead to an enhanced pressure. 
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It seems appropriate at this point to emphasize that SWa 7^ PSLp where as above, 6Lp = rjSL is the total length 
variation of the rods: expanding a given rod differs from moving a piston in that the piston has to be held fixed, while 
our move pertains to a given particle in the fluid, that is therefore free to move. If we would consider the particle 
expanded to be fixed at a given location, and with a slowly increasing size {£ —^ £ -\- Si), the work required would 
reduce to PSi since the particle would act as a piston. This work can again be expressed as the sum of a kinetic 
pressure term with the same value of the pair distribution at contact as given by Eq. I113\) . and an energy potential 
difference: 

PS£ = Pcont.ctkTSe + 5{Njr*p) (21a) 

= p kTSl + 5(Njr* p). (21b) 

1 — 77 

The difference with Eq. p9)) is that now the range r* is fixed but the density varies according to Jp/p ~ —5L/L = 6£/L. 
Hence, 

PS£ = —^kTSi + jr* p^S£, (22) 
1-ri 

and we recover Eq. In absence of the Kac tail (7 = 0), fixing or not the expanded particle is unimportant and 

the fact that we there have SWa = PSLp can be viewed as a consequence of the hard potential used (with interactions 
at contact only and a pressure that only depends -apart from p and T~ on the pair distribution function at contact) 



IV. CONCLUDING REMARKS 



We have proposed an exact derivation of the equation of state of certain simple liquids: hard rod systems with 
and without a long range pair potential of interaction (Kac tail). The argument -based on scaling considerations- 
can be easily extended to hard core particles in higher dimensions (hard discs, hard spheres. . . ). The interest of the 
approach lies in its simplicity and it is our hope that the method put forward here is instructive and may be useful 
to illustrate an advanced undergraduate statistical mechanics course. 

To our knowledge, pegagogical accounts avoiding mathematical complication on the present topics are scarce in the 
literature. We note however that a relatively simple approach generalizing the original Bernouilli derivation of the 
ideal gas equation of state has been proposed, [isj In addition to being heuristic, this method could not be transposed 
to the case of Kac tails. 

Finally, we note that our arguments bear some similarities with the scaled particle theory developed in the 1960s, 
tiut are nevertheless distinct. Scaled particle theory is based on the computation of the reversible work 
needed to create a spherical cavity from which the centers of other spheres are excluded. This expression is then 
related to the density of particles at contact with the cavity boundary. This can be achieved exactly in one dimension, 
since there are then no curvature effects. The latter remark also explains why we have been able to derive the pressure 
in closed form for d — 1. 
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